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PROCEDE DE C&RACTERI SAT ION OPTIQUE DE MATERIAUX SAMS 
UTILISATION DE MODELE PHYSIQUE 

DESCRIPTION 

5 DOMAINE TECHNIQUE 

La presente invention concerne un procede 
de caracterisation optique de materiaux. 

Ce procede permet de caracteriser des 
couches minces ou epaisses de ces materiaux, qui sont 
10 formees sur des substrats. Les grandeurs physiques, que 
ce procede permet de determiner, sont: 

- l'epaisseur d'une couche d ! un materiau, 

- l'indice de refraction de ce materiau, et 

- le coefficient d 1 absorption de ce 

15 materiau- 

La caracterisation optique des materiaux 
est utile pour 1 T analyse chimique de ces materiaux 
(notamment etude des bandes d* absorption, des 
proprietes de densif ication et des proprietes 

20 d' oxydation) , dans les domaines de la micro- 
electronique, des capteurs, de la biologie, de la 
medecine) , ou bien pour l 1 analyse des epaisseurs de 
depots de ces materiaux . 

Pour des exeraples d ! application, on se 

25 reportera au document [1] qui, comme les autres 
documents cites par la suite, est mentionne a la fin de 
la presente description* 

La cciracterisation des proprleles optiquos 
d'un materiau est egalement utile lorsque le materiau 

30 est par la suite structure ( pour y former par example 
des gravures ou des rugosites) et que les proprietes de 
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diffraction optique de la structure obtenue doivent 
etre calculees (voir le document [2]). 

Indiquons des a present que 1 T invention est 
particulierement utile lorsque la loi physique suivie 
5 par l'indice de refraction complexe du materiau que 
l'on veut caracteriser n T est, a priori , pas connue- 

ETAT DE IA TECHNIQUE ANTERIEURE 

On rappelle que les mesures optiques 
10 peuvent etre de diverses natures : 

II peut s'agir de mesures 

ref lectometriques . Dans ce cas, le coefficient de 
reflexion en intensite d'une structure est mesure sur 
un spectre (c 1 est-a-dire un intervalle) de longueurs 
15 d'ondes [Xn, X M ] * 

L T angle d r incidence de la lumiere 
d T eclairement peut etre non nul* Le coefficient de 
reflexion peut etre mesure pour plusieurs angles 
d ! incidence 0. Nous noterons R (6, X f p) le spectre 
20 ref lectometrique, ou p est la polarisation du faisceau 
incident et X la longueur d'onde de ce dernier. 

Generalement, 1' angle 9 est nul et la 
polarisation p indeterminee . Dans le cas ou 9 n'est pas 
nul, il faut connaitre cette polarisation p. En 
25 general, cette derniere est de type (S) ou (P) . 

II peut s'agir aussi de mesures 
ellipsometriques . Les grandeurs mesurees sont alors les 
parties reelle et imaginaire du rapport du coefficient 
de reflexion en polarisation (P) au coefficient de 
30 reflexion en polarisation (S) • 
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On note genera lement p=lp|exp(jA) ce 
rapport complexe (avec j 2 =-l) . Les grandeurs 
generalement traitees sont |p|, que l'on note tan(\|/), 
et cos (A) f ou des combinaisons des deux. 
5 Par exemple, les variables issues d'un 

ellipsometre a modulation de phase sont les suivantes : 
I s =sin (2\j/) sin (A) et I c =cos (2\[f ) . 
Un ellipsometre classique fournit, quant a 
lui, les variables suivantes : 
10 a=(tan 2 \j/-l) / (tan 2 \|/+l) et (3=cosA (1-a 2 ) 1/2 . 

Par souci de generality , nous noterons Si et 
S2 les grandeurs traitees. 

Les spectres Si, ie [1,2], sont mesures sur 
une plage de longueurs d T ondes [A™, X M ] * L f angle 
15 d ? incidence peut etre quelconque. Plusieurs spectres 
p euven t etre mesures a differents angles d' incidence 
afin d'obtenir un spectre plus riche. Nous noterons 
s(9,A)={Si (9, A,) f S 2 (Q f X)} le spectre ellipsometrique . 

De fagon complementaire, des mesures 
20 goniometriques (coefficient de reflexion en fonction de 
1' angle d' incidence) peuvent etre ajoutees aux mesures 
servant a la caracterisation, afin de determiner 
l'epaisseur des diverses couches, pour une ou plusieurs 
longueurs d'ondes. Ces mesures ne sont pas suffisantes 
2E3 en elles-memes puisque 1 ' on veut determiner 1 1 indice de 
refraction complexe sur une gamme spectrale allant de 
A™ a Am . 

Afin de sirnplifier 1' expose, nous noterons 
W un ensemble de spectres ref lectometrique ou/et 
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ellipsometrique (s) (et eventuellement goniometriques 
pour quelques longueurs d'ondes). 

Sans perdre aucune generality, nous 
n 1 exposerons, dans la presente description, le mode 
5 d 1 utilisation des precedes de l T art anterieur et de la 
presente invention que dans le cas d'une seule couche 
mince d T un materiau, formee sur un' substrat connu. 

L'epaisseur de cette couche est notee e et 
1 T indice de refraction complexe du materiau a la 
10 longueur d' onde X est note n* (X). 

On rappelle a ce propos que la partie 
reelle ( respect ivement imaginaire) de cet indice de 
refraction complexe est notee n (X) (respectivement 
k{X)) et appelee "indice de refraction" (respectivement 
15 "coefficient d' extinction") . 

En outre on note Er (^¥ n) , T (2) ) une fonction 
d'erreur (par exemple l'ecart quadrat ique moyen) entre 
deux spectres W a) et T (2) * 

Par exemple, on pourra prendre, lorsque 
20 l'on dispose de spectres ellipsometriques sur plusieurs 
angles Q lf ie {l..n}, et d'un spectre ref lectometrique : 
E c C¥ m , X Y t2> ) = 

+ [r 0) (A) - R {2) (A)] 2 JrfA (1) 

25 

avec (X)={Si n) (6 £ ,A),S 2 (1) (B lf X) ,R (1) (X)} 

¥ <2) {X) ={Si {2) {Qi,X) , s 2 (2) (QuX) ,R {2) (X)} 
et ie{l...n} 
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Des facteurs de ponderation peuvent eLre 
apportes a l'integrale de fagon a ce que la fonction 
d'erreur puisse tenir compte de variations sur la 
precision de mesure des spectres. 
5 La caracterisation optique de couches de 

materiaux s'articule generalement autour de deux 
applications : 

La premiere application est le controle 
dimensionnel du depot de couches minces que l'on 

10 utilise en microelectronique . 

Generalement on connait bien le materiau 
depose, c f est-a-dire que 1 1 on connait bien 1 1 indice de 
refraction complexe de ce materiau aux longueurs 
d'ondes de la lumiere utilisee pour la caracterisation. 

15 Les lois suivies par l 1 indice de refraction 

complexe sont soit tabulees, soit approchees par des 
lois physiques connues telles que, par exemple, le 
modele de Cauchy, le modele de Sellrneier (voir le 
document [3] ) , les lois de Forouhi (voir le document 

20 [4]), et les lois lois d T oscillateurs harmoniques (voir 
le document [5]). Ces lois sont definies par un nombre 
fini de paramctres. 

Par exemple, une loi de type Cauchy sans 
absorption, a deux parametres, est definie de la 

25 maniere suivante : 

Refn*U)] =n(X)=a 0 + ^V 
Im[>i *(A)]«k(A.)«0 

Lorsque 1 T on est sur de la valeur des 
coefficients aide {0,1}) mais que l'on ne connait pas 
30 l'epaisseur, un algorithme de recherche est utilise 
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afin de trouver 1 1 epaisseur qui minimise l'erreur entre 
la mesure Y et la reponse theorique V compte tenu de 
l 1 indice modelise. 

L'algorithme de recherche peut etre, par 
5 exemple, la methode du Simplex, la recherche Tabou, la 
methode de Levendt-Marquart ou la methode du recuit 
simule (voir le chapitre 10 du document [6]). 

Lorsque 1 1 indice de refraction est 
approximatif , les coefficients a± sont integres dans la 
10 procedure d'ajustement de Y et La recherche des 

coefficients ax constitue un procede de caracterisation 
de l 1 indice de refraction . 

Cependant, lorsque la loi suivie par cet 
indice de refraction est inconnue (il arrive que le 
15 materiau soit inconnu ou qu'il ne soit pas bien decrit 
par une loi physique connue) , ce procede reste 
approximatif et l 1 epaisseur risque d'etre fausse. 

La deuxieme application est la 
caracterisation des materiaux. 
20 Le procede utilise reste le meme, mis a 

part que le materiau n'est pas bien connu. C'est 
precisement la fonction d' indice de refraction complexe 
la plus proche de la realite qui est visee. 

Le type de loi peut etre choisi par 
25 analogie avec d'autres materiaux. Cependant, la loi 
suivie par 1 T indice de refraction complexe peut etre 
compliquee, ce qui est par exemple le cas d'une loi 
d T oscillateurs harmoniques : 

[n(B) + jkQBtf = 1 + V * dL 

+ jTi E~Ei + jTi 
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Dans 1 1 expression ei-dessus, j "=-1 et 
l'indice de refraction et le coefficient d' extinction 
sont exprimes non pas en fonction de X mais de E, avec 
E=: 1PA0/X (X en nra) . 
5 Dans ce cas, les coefficients des 

osciliateurs sont difficiles a trouver si 1 ! on n ! a pas 
leur ordre de grandeur. La recherche est dif f icilement 
automatisable, les algorithmes de recherche pouvant 
donner des reponses erronees et le temps perdu peut 

10 etre considerable. 

II existe une alternative & la recherche 
des coefficients : la methode point~a-point (PAP) . 
Cette methode PAP propose de ne pas choisir de loi 
physique et de rechercher l'indice de refraction 

15 complexe du materiau pour chaque longueur d'onde X if ou 
i<= [l.«.n], avec Xi-X^ et A„^X M . 

Pour chaque Xj , un algorithme de recherche 
tente de trouver l'epaisseur, l'indice n {Xi) et le 
coefficient d T extinction k (X L ) qui minimisent l'erreur 

20 entre la mesure Y et la reponse theorique 

{X ir n{Xi) ,k(X k ) ,e) . 

Un tel procede pose un probleme parce que 
les divers points {X ir €, n (Xi> , k {X x ) ) ne sont pas 
forcement physiquement compatibles entre eux : par 
25 exemple, repaisseur trouvee peut varier en fonction de 
la longueur d'onde et la loi suivie par l'indice de 
refraction complexe, plus s implement appelee loi 
d'indice, peut presenter des discontinuites . 
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Ce procede est generalement valable 
seulement lorsque I'epaisseur est tres bien connue et 
que les mesures sont de tres bonne qualite. 

EXPOSE DE L' INVENTION 

La presente invention a pour but de 
reraedier aux inconvenients precedents . 

Le procede objet de 1' invention permet de 
caracteriser un materiau. sans utiliser un modele 
physique, c'est-a-dire sans utiliser une loi physique 
suivie par 1' indice de refraction complexe du materiau 
etudie. II est done tout particulierement utile 
lorsqu'une telle loi n'est pas connue. 

Ce procede constitue une alternative aux 
procedes de caracterisation connus, mentionnes plus 
haut. II peut etre appele "methode des noeuds" car il 
utilise des "noeuds" e'est-a-dire des points de 
coordonnees (Ai n ] ) r ou n* est la valeur prise par 

l f indice de refraction complexe a la longueur d' onde \± 
et 1 prend un nombre limite de valeurs (entieres) . 

De fagon precise, la presente invention a 
pour objet un procede de caracterisation optique d'au 
moins une couche d'un materiau dans un intervalle A de 
valeurs prises par une fonction a d'une longueur d' onde 
optique X, lorsque X vari'e dans un intervalle de 
longueurs d'ondes, cette couche §tant formee sur un 
substrat, ce proced6 etant caracterise en ce qu'il 
comprend les etapes suivantes : 

1) on effectue un ensemble de mesures de 
ref lectometrie et/ou d T ellipsometrie sur 1' intervalle 
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A, eel ensemble de mesures conduisanL a un spectre 
rnesure, note Y, et l'on choisit les methodes de calcul 
associees a la nature des mesures et au type de couche 
a caracteriser, 

5 2) on choisit m valeurs initiales a 3 ...o^ de 

la fonction a, appartenant a cet intervalle A, m etant 
un nombre entier au moins egal a 1, et l'on definit un 
intervalle B comme etant 1' ensemble des points a de 
1' intervalle allant du plus petit au plus grand des 

10 nombres ai„.Oni, lorsque m est superieur a 1, et comme 
etant 1' intervalle A lorsque m vaut 1, 

3) on choisit m valeurs initiales complexes 
d'un indice de refraction complexe n^n+jk aux m points 
a i7 i allant de 1 a m, 

15 4) lorsque m est different de 1, on choisit 

une loi d' interpolation qui permet de calculer l f indice 
de refraction n(a) du materiau sur 1' intervalle B, a 
partir des points (a if n^) , avec ni=n (ai) , i allant de 1 
a m, et J orsque m vaut 1, n{a) est pris egal au nombre 

20 ni(a!) sur tout 1' intervalle B, 

5) on choisit M parametres variables, M 
etant inferieur ou egal a 2m+l, 

6) on choisit une fonction d'erreur 

Er f y y ) qui caracterise l'ecart entre un spectre 

25 mesure X Y el un spectre experimental V F , 

7) a I'aide d'une fonction de minimisation 
do Kr CV , X V ) a M parametres, on effectuo 1 * sorie 
d' e tapes suivante : 
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a) a l'aide de la loi d' interpolation des 
(a lf rii) sur 1' intervalle B, on deduit n(oc), a 
appartenant a B, 

b) a l'aide de n(a) et de l'epaisseur s de 
5 la couche, et de methodes de calcul de spectres, on 

calcule un spectre theorique Y (n (a) , e) , 

c) on compare Y et W a l'aide de Er C¥, W) 
et, si Er(Y, T) est suf f isamment petite, c'est-a-dire 
inferieure a une valeur predefinie e, ou est minimale, 

10 on passe a l'etape e) , sinon on passe a l'etape d) , 

d) on fait varier les M parametres 
variables de fagon a tendre vers le minimum de Er OF, 
W ) et l'on retourne a l'etape a), 

e) si ErpF, Y) est inferieure a e, on 
15 obtient done un ensemble de M parametres variables, 

pour lequel Er C¥ , *F (n (a, M) , e) ) est minimum et l'indice 
de refraction est alors pris egal a celui qui a ete 
obtenu en dernier, et si Er C¥ f *F ) est superieure ou 
egale a e on passe a l'etape 8), 
20 8) on augmente le nombre m de valeurs 

initiales de la fonction a et l'on retourne a l'etape 
2) . 

II est done par exemple possible 
d'effectuer la caracterisation optique : 
25 - sur un intervalle de longueurs d' ondes X, 

en 1' occurrence 1' intervalle [Xmin, Xmax] , 

- ou sur un intervalle de longueurs d' ondes 
inverses 1/^, en 1' occurrence sur un intervalle 
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[(1/X)min, (1/A.)max] , ou (1/X)rnin est egal a 1/ (Xmax) 
et (lA)max a l/(Amin), 

ou sur un intervalle d' energies E (avec 
E = hv = hc/X ou h est la constante de Planck, c la 
5 vitesse de la lumiere dans le vide et v la frequence 
correspondant a X) , en 1' occurrence sur un intervalle 
[Emin, Emax] , ou Emin est egal a he/ (Xmax) et Emax a 
he/ (Xmin) , 

ou, plus generalement, sur un intervalle 
10 [amin, amax] de valeurs prises par une fonction a de la 
variable X. 

II convient en outre de noter que 
1' invention est utilisable pour caracteriser un spectre 
ou une partie de spectre, 

15 Chaque loi d 1 interpolation peut etre 

choisie parrni les lois d 1 interpolation lineaires, les 
lois d' interpolation cubiques, les lois d 1 interpolation 
polynomiales et les lois d f interpolation par exemple de 
type fonction « spline » . 

20 Selon un mode de mise en oeuvre prefere du 

procede objet de 1' invention, le spectre est 
echantillonne uniformement en a ( X) , e'est-a-dire que 
les valeurs initiales de la fonction a (voir l'etape 2) 
rnentionnee plus haut) sont reparties de maniere 

25 uniforme dans 1' intervalle A, la repartition des nceuds 
etant ainsi homogene , 

Comme on l 7 a vu f a (A.) peut etre choisi 
parmi X, 1/X et hc/X ou toute autre fonction de X, ou h 
est la constante de PlancI: et c est la vitesse de la 

30 lumiere dans le vide, 
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De preference, a 1'etape 6) mentionnee plus 
haut, on mesure 1'erreur sur un intervalle d' inter^t C 
qui est inclus dans 1' intervalle B ou egal a cet 
intervalle B. 

5 Les M parametres variables peuvent etre les 

parties reelles des indices de refraction aux points 
cti, i allant de 1 a m, ou les parties imaginaires de 
ces indices de refraction, ou ces M parametres 
variables peuvent etre constitues par l'6paisseur du 
10 materiau dont on cherche l'indice de refraction* 

BREVE DESCRIPTION DES DESSINS 

La presente invention sera mieux comprise a 
la lecture de la description d'exemples de realisation 
15 donnes ci-apr&s, titre purement indicatif et 

nullement lirnitatif, en faisant reference aux dessins 
annexes, sur lesquels : 

- la figure 1 est une vue schematique de 
dispositifs permettant de caracteriser une couche 

20 conformement a 1' invention, 

la figure 2 montre les variations de 
l'indice de refraction en fonction de la longueur 
d'onde, pour un materiau suivant une loi de Cauchy 
(courbe I) et pour un materiau caracterise conformement 

25 & 1' invention (courbe II), 

- la figure 3A (respectivement 3B) montre 
les variations de l'indice de refraction 
(respectivement du coefficient d' extinction) en 
fonction de la longueur d'onde, pour un materiau 

30 suivant une loi a deux oscillateurs harmoniques (courbe 
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I) et pour un rnateriau caracterise con! ormernent a 
1' invention, et 

- la figure 4 illustre schernatiquement les 
param^tres utilises dans une generalisation d'exemples 

5 de 1' invention. 

EXPOSE DETAILLE DE MODES DE REALISATION PARTICULIERS 

L r invention propose une alternative aux 
procedes classiques, mentionnes plus haut . Elle permet 

10 d'allier la coherence d'un modele de couche 
(correspondant a une loi d'indice continue et a des 
epaisseurs physiques constantes) , a la generality 
concernant la loi d'indice a trouver (comine dans la 
methode PAP). En outre, la resolution n'est limitee que 

15 par la resolution du spectre mesure. 

Dans le procede objet de 1' invention, le 
spectre d'indice n* (A) est caracterise par : 

- un nombre reduit de "nceuds", qui sont des 
points de coordonnees (Aj,nj,ki) ou (Aj, n*j), avec ni = n 

20 (Ax), ki « k (A x ) et n'j^ni+jki, ou j 2 =-l, et 

- une loi d T interpolation entre les nceuds, 
qui peut etre, par exemple, lineaire, cubique, de type 
« spline » ou polynomiale (de degre quelconque) . 

Cette loi d 1 interpolation permet de 
25 calculer, 3 partir des noeuds, les indices de refraction 
et les coefficients d 1 extinction pour les longueurs 
d 1 Grides situees entre les noeuds- 

Par exemple, lorsque 1 1 indice de refraction 
est caracterise par un ensemble de valeuus pour des 
30 longueurs d'ondes A) ... A m , on peut utiliser une 
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interpolation lineaire entre deux longueurs d'ondes \± 
et Xi+i pour calculer 1 ' indice n a la longueur d'onde X 
(voir le document [6] chapitre 3) : 



On peut faire de meme pour le coefficient 
d' extinction . 

Lorsque le nombre de noeuds le permet, des 
formules d 1 interpolation plus complexes, faisant 
10 intervenir les nceuds voisins, peuvent etre utilisees 
(voir le document [6] chapitre 3) . 

Un raodele de couche est done caracterise 
par une epaisseur e et une famille de noeuds. 

On donne ci-apres un exemple du precede 
15 objet de I 1 invention. Dans cet exemple, les mesures Y 
sont constitutes d'une mesure ref lectometrique R (X) et 
d'une mesure ellipsometrique Si, 2 (©,^0 ou B est l'angle 
d f incidence du faisceau lumineux que l'on envoie sur la 
couche a etudier lors de la mesure d 1 ellipsometrie . 
20 Cette couche est une couche mince de sorte 

que 1' tpaisseur de cette couche est aussi une variable 
du probleme. En outre, on suppose qu'une seule couche 
est inconnue, cette couche etant formee sur un substrat 
connu. 

25 Expliquons d'abord brievement cet exemple, 

qui utilise un algorithme (algorithme de la « methode 
des noeuds » conformement a 1' invention) . 

A partir d 1 informations supposees sur 
1' epaisseur e de la couche etudiee et sur l f indice de 



n(X) 




5 



avec \i<\<\i+i 



B 14382.3 PV 



15 

refraction n (X) et le coefficient d 1 extinction k {X) du 
materiau de cette couche, on construit des noeuds de 
depart (en faible nombre) et une epaisseur de depart E. 

On dispose ainsi de m noeuds et, par 
interpolation, on peut connaitre n (X) et k(X) en dehors 
des valeurs des longueurs d'ondes associees aux noeuds. 

A partir de 1' epaisseur de depart e et de 
ces valeurs de depart n (X) et k(X), on determine le 
spectre theorique W en utilisant des calculs 
ellipsometriques et ref lectrometriques . 

Par ailleurs, au moyen de dispositifs 
d'ellipsometrie et de ref lectrometrie et d ! un 
spectrometry, on obtient Si, a (&rX) et R(X) et l'on en 
deduit les mesures notees T (pour des conditions de 
mesure 0 et X) . 

On compare ensuite X F et \F en utilisant une 
fonction d'erreur Er et 1 1 on optimise la valeur de 
l'indice de refraction et la valeur du coefficient 
d 1 extinction aux differents noeuds, ainsi que la valour 
de 1* epaisseur, en cherchant a minimiser Er 

Lorsque ces valeurs sont optimisees et si 
la precision sur le spectre n(X) r le spectre k(X) et 
l'epaisseur e n'est pas suffisante, on ajoute de 
nouveaux noeuds, on fait varier I 1 epaisseur 8, et 1 1 on 
recommence la determination de \F , la comparaison de *F 
et W et 1 1 optimisation que l'on a mentionndes plus 
haut, etc. 

On arrete la boucle ainsi definie lorsque 
la precision sur chacun des spectres n [X) et k(A) et 
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sur l'epaisseur e est jugee suffisante (ajustement 
satisfaisant de Y et W) . 

Les spectres n (X) , k{X) et l'epaisseur e 
sont ainsi caracterises. 
5 Dans un exemple donne a titre purement 

indicatif et nullement limitatif, on trouve une 
epaisseur e egale a 212, 3nm pour la couche. 

La figure 1 raontre de fagon schematique la 
couche etudiee 2, formee sur un substrat 4. On voit le 

10 dispositif d' ellipsometrie 5,6, le dispositif de 
r6f lectrometrie 8 et le spectrometre 10. On voit en 
outre des moyens electroniques de traitement 12, 
comprenant un ordinateur et permettant de caracteriser 
n(X) , k{X) et 8 en fonction des informations fournies 

15 par le spectrometre 10 et conf ormement au proc^de de 
1 1 invention . 

Ces moyens 12 sont munis de moyens 
d'affichage 14 qui permettent, en particulier, 
d'afficher la courbe des variations de n en fonction de 
20 X et la courbe des variations de k en fonction de X. 

Revenons de fa<?on plus detaill«§e sur 
1 1 exemple donn6 . 

Phase 1 

Le proc6de de cet exemple comporte d'abord 
25 une- etape d' initialisation. 

L'algorithme debute avec un nombre reduit 
de nceuds, plus precisement au. moins un nceud. On peut 
done d6buter par un seul noeud, en imposant un indice de 
refraction et un coefficient d 1 extinction qui restent 
30 constants lorsque la longueur d'onde varie. 
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On choisit des positions cle nceuds cie fa^on 
a pouvoir, a partir de cette famille de nceuds, deduire 
tout le spectre par interpolation. Le module de couche 
est done a 3 parametres ou plus, puisque l'epaisseur 

5 est aussi une variable a determiner. La table d'indice 
sur tout le spectre est done deduite des nceuds par 
interpolation . 

C'est le cas lorsque l'epaisseur des 
couches est de 1 1 ordre de grandeur de la longueur 

10 d'onde (couches minces). Mais lorsque l'epaisseur de la 
couche est superieure a la longueur de coherence de la 
source de lumiere, l'epaisseur n'intervient quasiment 
plus dans le calcul de la reponse de la couche, et 
n'est done plus une variable du probleme. 

15 Par exemple, dans le cas de disques 

optiques (CDROM) , sur lesquels une couche tres 
importante est deposee (son epaisseur est de 1 ! ordre de 
1 millimetre), le coefficient de reflexion d'une telle 
couche n f est plus fonction de l'epaisseur de la couche 

20 mais seulement de l'indice de refraction de celle-ci, 
la longueur de coherence du faisceau lumineux incident 
etant inferieure a l'epaisseur de la couche. Dans ce 
cas precis, la longueur de coherence du faisceau 
incident est determinee par la rugosite des couches. 

25 On choisit, par exemple, de placer les deux 

premiers nceuds aux extremites X min et X max du 
spectre. Les valeurs de 1 1 indice complete en ces 
extremites sont choisies en fonction du type du 
materiau etudie. Par exemple, sur un spectre 

30 ellipsornetrique entre 300nm et 300nm d'une couche mince 
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de resine photosensible, on prend n (300nm) =n (800nm) =1, 5 
et k(300nm)=k(800nm)=0. 

Lorsque le spectre n'est caracterise que 
par deux noeuds, l'indice entre les extremes est 
5 determine par interpolation lineaire. Dans le cas 
considere, on a done n(X)=l,5 et k(X)=0 pour X 
appartenant a [300nm, 800nm] . 

A partir de trois noeuds, on choisit plutot 
une interpolation cubique afin d'obtenir des formes de 
10 loi d f indice plus douces que celles obtenues par 
interpolation lineaire . 

L f epaisseur de depart est, quant a elle, 
choisie aussi proche que possible de l'epaisseur. 
reelle . 
15 Phase 2 

On precede ensuite a une determination 
optimale des valeurs de l T indice de refraction et du 
coefficient d T extinction sur les noeuds et de la valeur 
de 1 T epaisseur. 

20 Pour ce faire, les spectres Y (A) sont 

calcules a l'aide " du modele de couche utilise, 
resultant du choix des noeuds, de la loi d 1 interpolation 
et de 1' epaisseur de la couche. 

Le modele physique utilise pour le calcul 

25 de W est bien sur fonction de la methode de mesure 
utilisee, c r est-a-dire notamment de' 1 1 angle d' incidence 
de la lumiere, du spectre utilise et du modele. de 
couches minces ou de couches epaisses le cas 6cheant 
(voir par exemple les modeles de couches empilees dans 

30 le document [3] ) . 
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Les spectres H' etant constitues par un 
ensemble de mesures de natures diverses (par exemple 
des mesures ellipsometriques et ref lectometriques) , on 
utilise pour les mesures de ref lectometrie 
5 (respectivement d' ellipsometrie) une methode de calcul 
de ref lectometrie (respectivement d 1 ellipsometrie) • 

Les mesures ref lectometriques et 

ellipsometriques sont combinees par 1 ' intermediaire 
d'une fonction d'erreur Er('SK / x F) qui est par exemple du 
10 genre de celle qui est definie par I'equation (1). 

Grace • a une fonction de recherche , on 
minimise l'ecart entre W (X) et T (X) , en faisant varier 
la valeur de I'indice de refraction et la valeur du 
coefficient d f extinction a la position de chacun des 
15 nosuds ainsi que 1 ! epaisseur de couche (si cette 
epaisseur est un facteur important dans le calcul de 
W ) . Lorsque l'ecart est minimum c f est-a-dire lorsque 
Er( x F,Y) est minimum, cela signifie que les mesures 
ref lectometriques et ellipsometriques coincident au 
20 mieux (pour un nombre de noeuds donne) . 

A ce stade, on obtient, pour un nombre 
connu de noeuds et une position spec t rale connue pour 
chacun de ces des noeuds, le model e de couche (indice de 
refraction, coefficient d 1 extinction et epaisseur) qui 
25 correspond au rnieux a la couche reelle. 

La validite du modele trouve est d'autant 
plus assuree que le nombre de mesures est grand. Pour 
avoir un grand nombre de mosures, on pout par oiemplo 
utiliser plusieurs angles d T incidence de la lumiere 0.i, 
30 1 i i ^ *' f pour les spectres ellipsonielciques, fa ire 
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une mesure ref lectometrique et faire des mesures 
goniometriques complementaires . 
Phase 3 

Ensuite, on augmente le nombre de nceuds. 
5 On ajoute un nombre fini de nceuds. Dans un 

premier mode de realisation , les nceuds rajoutes sont 
positionnes de fa<?on & representer au mieux les 
spectres n{\) et k(X). A titre d'exemple, on place ces 
nceuds supplementaires aux endroits ou l T ecart entre Y 

10 et W est maximum ou aux endroits en lesquels les nceuds 
sont les plus espaces. Et 1 1 on retourne a la phase 2 
tant que la precision sur chacun des spectres n (X) et 
k(X) et sur l'epaisseur e n'est pas suffisante r 
c'est-a-dire n'est pas egale a une precision 

15 predefinie. 

Dans un deuxieme mode de realisation, on 
intercalle de nouveaux nosuds entre deux noeuds de 
1' ensemble de nceuds prealablement choisis, ces nouveaux 
noeuds se repartissant de fa?on homogene sur le 
20 spectre. 

II convient toutefois de noter ce qui suit. 
Lorsque l'on augmente le nombre de noeuds, 
il n'est pas necessaire de conserver la position des 
anciens noeuds. Par exemple, admettons que l'on 

25 echantillonne uniformement un spectre de 400nm a 800nm 
avec 3 nceuds. Ces noeuds se situent done respectivement 
a 400nm, 600nm et 800nm. Lorsque l'on passe a 6 nceuds, 
les trois noeuds additionnels ne peuvent etre places de 
fa<jon a ce que le spectre soit uniformement 

30 echantillonne si l'on conserve la position des anciens 
noeuds. La position des 6 noeuds peut etre d6finie, si 



B 14382.3 PV 



21 



l'on souhaite un echantillonnage uniforme par les 
valeurs 400, 480, 560, 640, 720, 800nm. L'ancien noeud 
central a 600nm clisparait done. Pour calculer la valeur 
de l'indice a ces positions a partir des anciens nceuds, 
5 on procede par interpolation. 

Dans ce qui suit, nous donnons deux 
exemples courants d 1 application de l 1 invention. Ces 
exemples mettent en jeu deux types de materiaux qui 
suivent des lois differentes. 

10 A partir de mesures ellipsometriques et 

ref lectometriques, nous proposons de retrouver les lois 
physiques suivies par ces materiaux. 

Nous procedons de la fa?on suivante : 

Un materiau fictif est cree, ce materiau 

15 suivant une loi theorique connue (une loi de Cauchy ou 
une loi d 1 oscillateurs harmoniques) , avec des 
parametres que nous fixons arbit rairement . Les 
variations de l ! indice de refraction complexe en 
fonction de la longueur d'onde sont ainsi parfaitement 

20 connues . 

Nous imposons en outre au materiau unc 
epaisseur de 200,00nm, sur un sufostrat de silicium, ce 
dernier etant egalement tres bien connu. 

Des mesures fictives (mesures 

25 ellipsometriques, mesures ref lectometriques) sont 
calculees, puis bruitees de fagon a introduire un 
def aut d 1 appareillage . 

Tout se passe conime si nous avions des 
mesures replies, effectuees sur le materiau. Mais, 
30 contrairement a la realite, nous connaissons 
parfaitement 1 1 indice de reftaclion complexe puisque 
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nous l'avons fixe, de meme- que nous avons fixe 
1' epaisseur de la couche de materiau. 

Nous testons ici la methode "en aveugle", 
c r est-a-dire que nous partons d'une epaisseur fausse 
(220nm) et d f indices de refraction complexes faux, 
puisqu'ils sont senses etre inconnus. 

Nous appliquons le procede objet de 
1' invention puis nous comparons l'indice de refraction 
complexe trouve a I'indice de refraction complexe 
theorique. Nous retrouvons bien les memes lois, tres 
precisement, ainsi que la meme epaisseur de couche. 

Comme premier exemple, prenons un materiau 
dont I'indice de refraction complexe suit une loi de 
Cauchy telle que : 

300 2 300 4 
»(A) = 1,5 + 04^ + 04^ 

k(X)=0 

Cette loi d'indice est typique des resines 
photosensibles (dans la plage spectrale allant de 300nm 
a 800nm) . 

Afin de retrouver a 1 1 aide de la methode 
des nceuds (c 1 est-a-dire du procede objet de 
l 1 invention) la loi d T indice mentionnee ci-dessus, nous 
effectuons deux mesures, a savoir une mesure 
ellipsometrique a un angle de 70° et une mesure 
ref lectometrique ♦ 

Les conditions de traitement sont les 

suivantes : 

- le spectre traite est compris entre 300nm 

et 800nm, 
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- au depart, les nceuds sent au>; positions 
(300nm, 1,6) et (800nm, 1,6), e'est-a-dire que l'indice 
est considere comme variant lineairement entre 300nm et 
800nm, et que sa valeur est constante (egale a 1,6), 

5 le nombre de nceuds est augmente 

iterativement selon la sequence 2— »4— >6, 

- lors de la procedure d' augmentation des 
nceuds, la position en longueur d' onde de chaque nceud 
est calculee de fa<?on que 1' echantillonage en 1/X soit 

10 uniforme (X : longueur d'onde), la densite -de points 
<§tant done d'augmentee vers les faibles longueurs 
d' ondes, 

- la loi d 1 interpolation est une loi 
cubique, lorsque le nombre de nosuds est superieur a 2, 

15 sinon elle est lineaire, et 

- l'algorithme de minimisation utilise est 
un algorithme de type Simplex. 

La figure 2 permet de comparer l'indice de 
refraction correspondant au materiau fictif qui suit 
20 parfaitement uno loi de Cauchy (courbe I) a l'indice de 
refraction que nous trouvons par la methode des nceuds 
(courbe II), a 1 1 aide de 6 nceuds (representees par des 
cercles sur la figure 2) . Nous avons utilise un X Y 
compose d'une mesure ellipsometrique {Si (X) , S 2 (X) } a 70° 
25 et d'une mesure ref lectometrique R(X). 

Nous about issons a une £paisseur de 

199, 8384nm. 

Comme deuxieme oxemple, prenons un matcr.iau 
dont l'indice de refraction complete suit une loi a 
30 deux oscillatenrs harmoniques, telle que : 
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avec 

j 2 = -1 et E = 1240/A (X en nm) 
Ai=0 , 25x1 , 5 2 A 2 =0 r 25x1, 5 2 
5 E L =1240/400 E 2 =12 40/300 

Gi=0,3 G 2 =0,3 

Dans ce deuxieme exemple, la methode des 
noeuds est appliquee a un ensemble de mesures 
ellipsometriques effectuees entre 250nm et 800nm, a 

10 75°, 70° , 60° et 45°, avec en plus une mesure 
reflectometrique. L'epaisseur reelle du materiau etant 
de 200nm, on trouve une 6paisseur de 200,25 nm avec la 
methode des noeuds, L'ajustement sur la loi d 1 indice 
eonsideree dans ce deuxieme exemple est tres bon, comme 

15 le montrent les figures 3A et 3B . 

Ces figures 3A et 3B illustrent 
respectivement les reconstructions des courbes n (X) et * 
k(X) du materiau a l T aide de la methode des noeuds. La- 
reconstruction est effectuee 1 1 aide de quatre 

20 spectres ellipsometriques et d'un spectre 

reflectometrique. Les pics d 1 absorption reels sont tres 
bien representees par la courbe obtenue par 
interpolation cubique entre les ncexids (representees par 
des cercles sur les figures 3A et 3B) . 

25 Sur la . figure 3A, la courbe I 

(respectivement II) correspond a un indice de 
refraction n qui suit parfaitement la loi choisie 
(respectivement a un indice de refraction n trouve par 
la methode des noeuds) . 
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Sur la figure 3B, la courbe I 
(respectivement II) correspond a un coefficient 
d' extinction k qui suit parf aitement la loi choisie 
(respectivement a un coefficient d' extinction k trouve 
5 par la methode des noeuds) . 

On vient de decrire des exemples de 
1' invent ion . On remarquera d'une maniere plus generale 
que, dans cette derniere, on considere un ensemble X de 
valeurs, avec X= {n 2 r n 2 , n ie ...n m/ k 3 r k?, ki, k m/ e} , ou 
10 n A est la valeur de I'indice de refraction 

(reel) au noeud correspondant a X if ie{L.m} r m etant le 
nombre de noeuds 

ki est la valeur du coefficient d 1 absorption 

au nceud correspondant a Xi, ie{l...rn}, 
15 e est l'epaisseur de la couche etudiee. 

Dans ce cas r l 1 operation de minimisation 
d'erreur Er C¥,W) revient a trouver 1' ensemble ou 
« vecteur » X tel que Er soit minimum. 

Lorsque 1 1 on n' impose pas de contrainte 
20 particuliere, la minimisation est une minimisation a 
2xm+l parametres. On peut bien sur ajouter des 
contraintes afln de diminuer le nombre de variables. 

En particulier, si 1' on sait que le 
materiau est non absorbant, on impose ki=0 pour tout i 
25 de {l...m} et X devient : X={"i , n*..., nj , ...n m , e} . 

Si, par une mesure complement aire (par 
exemple uno mesure de goniometrie ou une mesure directe 
non optiqu*) , on connait 1 1 ipaisseur de la couche 
consideree avec une precision suffisante, l f epaisseur £ 
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n'est .plus une variable et 1 1 on a : 
X={ni, Ti2r ~. f ni, ...n m , ki, k2, ki, »-r k m } . 

Bien entendu, les deux options precedentes 
peuvent etre combinees . 

On explique dans ce qui suit un mode de 
mise en ceuvre de 1' invention sous une forme plus 
generale que les exemples donnes precedemment . 

Soit a (X) une fonction de la longueur 
d'onde X de la lumiere utilisee pour les mesures. On 
peut par exemple choisir : 

a(X)=X (voir les figures 3A et 3B ou le 
spectre est unif ormement echantillonne en X) 

a(X)=l/X (voir la figure 2 ou le spectre 
est uniformement echantillonne en 1/X) 

a(X)=hc/X ou h est la constante de Planck 
et c la vitesse de la lumiere dans le vide, a(X) etant ■ 
alors homogene a une energie . 

Soit A l'intervalle du spectre de mesure, B 
1'intervalle spectral decrit par les « noeuds » et C 
l'intervalle d'interet. 

L'intervalle C est inclus dans l'intervalle 
B ou egal a cet intervalle B. De meme, l'intervalle B 
est inclus dans l'intervalle A ou egal a cet intervalle 
A. 

On precise que chacun des intervalles A, B 
et C est du type [ow o^j ou est inferieur a a M et il 
existe deux longueurs d' ondes X k et Aa telles que 
a m =a(X k ) et a^=a(X x ) • A titre purement indicatif et 
nullement limitatif, la figure 4 montre un exemple des 
intervalles A, B et C et des courbes des variations de 
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Y et de X F en fonction de a(X), Y etant en fait une 
fonction de n*(a(X)). 

Les cercles N representent les nceuds • 
n* est l'indice complexe, qui est exprime 
5 ici en fonction de cx(X) et dont les parties reelle et 
imaginaire sont respectivement notees n(a(A)) et k 

(a an . 

On voit egalement, a bitre purement 
indicatif et nullement limitatif, un exemple de la 
10 courbe des variations de n (respectivement k) en 
fonction de a (Pi) , passant par des points de coordonnees 
(n ir o&i) (respectivement (ki, a t ) , ou ai=a(Aj.), ni=n(ai), 
ki=k(ai), l<i<m (m entier naturel non nul) . 

Sur la figure 4, on notera la 
15 correspondance entre les noeuds N, les points (n ir a*) et 
les points (kj, ax). 

Dans le mode de mise en ceuvre considere, on 
utilise un algorithme comprenant les etapes suivantes : 

1. on fait les mesures sur 1' intervalle A 
20 et l'on choisit les methodes de calcul associees aux 

mesures (calculs ellipsometriques ou 

ref lectometriques) ; 

2. on choisit m n ombres oti (constituant m 
valeurs initiales de la fonction a) , i appartenant a 

25 {l,...,m}, avec m>l , et {o^} c A (les a t correspondant aux 
« ncouds ») ; lorsque m>l, B est defini cornme 1' ensemble 
des points a tels que min(ai) <a< rnax(Oi) / lorsque nv=l r 
on a B^A ; 
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3. on choisit m valeurs initiales d' indice 
complexe n* aux m points Ol, i appartenant a {l,...,m} ; 

4. si m^l on choisit une loi 
d' interpolation qui permet de calculer 1' indice de 
refraction n(a) sur l'intervalle B a partir des points 
(a*, ni), i appartenant a {l,,..,m} ; si m=l, alors 
n(a)=ni(ai) sur tout l'intervalle B ; 

5. on choisit M parametres variables avec 
M<2m+1 ; ces parametres peuvent etre par exemple : 

les parties reelles des indices de 
refraction aux points (Xi,i appartenant a {l,.-.,m}, ou 

- les parties imaginaires des indices de 
refraction en ces m§mes points, ou 

- l'epaisseur du materiau dont on cherche 
1' indice de refraction ; 

6. on choisit une fonction d'erreur 

Er OF, W) qui caracterise l'ecart entre un spectre 
experimental et un spectre mesure ; en general, k 
1'erreur est mesuree sur l'intervalle C ; 

7. a l'aide d' une fonction de minimisation 
de ErOF,^ 7 ) a M parametres, on effectue la serie 
d' etapes suivantes : 

a) a l'aide de la fonction d' interpolation 
des (CGi, ni) sur B, on deduit n(a) avec a appartenant a 
B ; 

b) a l'aide de n(a) de l'epaisseur e et des 
methodes de calcul de spectre, on calcule le spectre 
theorique W (n(a),e) ; 
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c) on compare Y et *F a l'aide de 
Er=Er (*¥ , W) ; si Er est suf f isamment petite (c'est-a- 
dire si Er est inferieure a une valeur predef inie e) , 
ou bien si Er est minimale, on passe a 1'etape e) , 

5 sinon on passe a 1'etape d) ; 

d) on fait varier les M parametres de fagon 
a tendre vers le minimum de Er et l'on revient a 
1'etape a) ; 

e) si Er est inferieure a e on obtient done 
10 un ensemble de M parametres tel que Er(Y,T (n(oc,M),e)) 

est minimum et le calcul de l'indice de refraction est 
acheve : cet indice est pris egal a celui qui a et6 
obtenu en dernier / et si Er est superieure ou egale a 
e, on passe a 1'etape 8) ; 
15 3. on augmente m et l'on revient a 1'etape 

2) . 

II convient de noter que la presente 
invention peut utiliser non seulemenb des 
echantillonnages en longueur d' onde (X) mais encore des 
20 echantillonnages en frequence (c/A, ou, plus simplement 
, en energie (hc/X) et, plus generalement, en 
parametre fonction de la longueur d'onde. 

II convient en outre de noter qu'une etape 
essentielle de 1' invention (etape 8) mentionnee plus 
25 haut) n'est pas limitee a l'ajout d'un ensemble de 
nceuds aux noeuds deja existants : elle englobe le cas 
plus general oil le nombre de nceuds augmente . Cela veut 
dire, dans un mode de realisation particulier de 
1' invention, qu'apres minimisation avec 3 nomcls, si 
30 l'on v«ut passeL a G nceuds au total, la position des 3 
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anciens nceuds est « effacee » de fa<?on, par exemple, a 
avoir une densite de noeuds constante sur le spectre* En 
pratique, c'est la meilleure option. 1/ information sur 
l'ancienne position des anciens noeuds n'est pas perdue 
car la valeur des indices aux anciens nceuds est 
utilisee pour calculer la valeurs des 6 nouveaux noeuds 
(en fait 3 nouveaux noeuds plus trois anciens noeuds) . 

Ainsi, selon un mode de raise en oeuvre 
particulier du procede objet de 1' invention, on peut 
augmenter le nombre de valeurs initiales de la fonction 
a en ajoutant une ou plusieurs valeurs aux valeurs 
initiales existantes ; mais, selon un mode de mise en 
ceuvre prefere, on peut augmenter le nombre de valeurs 
initiales de la fonction a en remplagant les valeurs 
initiales existantes par de nouvelles valeurs initiales 
dont le nombre est superieur au nombre des valeurs 
initiales existantes. 

La presente invention n'est pas limitee a 
la caracterisation de couches minces. Elle s' applique 
aussi a la caracterisation de couches epaisses. 

En outre/ la presente invention n'est pas 
limitee a la caracterisation d'une seule couche, formee 
sur un substrat. Elle s' applique aussi a la 
caracterisation de deux, ou plus de deux, couches 
formees sur un substrat. 

Les documents cit6s dans la presente 
description sont les suivants : 

[1] R.M.A. Azzam and N.M. Bashara, 
Ellipsometry and Polarized Light, North-Holland Physics 
Publishing, 1997, chapitre 6. 



B 14382.3 PV 



31 



[2] B. K . Minhas, S.A. Coulombe, S. Sohail, 
H . Naqvi and J.R. McNeil, Ellipsometric scatterometry 
for the metrology of sub-0 . 1-jam-linewidth structures, 
Applied Optics, 37(22) : 5112-5115, 1998. 
5 [3] M. Born and E. Wolf, Principle of 

Optics, Cambridge University Press edition. 

[4] A.R. Forouhi and I. Bloomer, Optical 
dispersion relations for amorphous semiconductors and 
amorphous dielectrics, Physical Review B, 34(10) : 
10 7018-7026, November 1986. 

[5] F. L. Terry, Jr., A modified harmonic 
oscillator approximation scheme for the dielectric 
constants of Al x Gai- x As, Journal of Applied Physics, 
70(1), 1991, pages 409-417. 
15 [6] W. H. Press, S. A. Teukolsky, W. T. 

Vetterling and B. P. Flannery, Numerical Recipes in C, 
Cambridge University Press, 1992, chapitres 3 et 10. 
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RE VEND I CAT I ON S 

1. Procede de caracterisation optique d'au 
moins une couche d ! un mater iau dans un intervalle A de 
valeurs prises par une fonction a d'une longueur d' onde 
5 optique A, lorsque X varie dans un intervalle de 
longueurs d'ondes , cette couche etant formee sur un 
substrat, ce procede etant caracterise en ce qu'il 
comprend les etapes suivantes : 

1) on effectue un ensemble de mesures de ref lectometrie 
10 et/ou d'ellipsometrie sur 1' intervalle A, cet ensemble 

de mesures conduisant a un spectre mesure, note *F, et 
l'on choisit les methodes de calcul associees a la 
nature des mesures et au type de couche a caracteriser, 

2) on choisit m valeurs initiales ai ...0^ de la fonction 
15 a, appartenant a cet intervalle A , m etant un nombre 

entier au moins egal a 1, et l'on definit un 
intervalle B comme etant 1' ensemble des points a de 
1' intervalle allant du plus petit au plus grand des 
nombres (Xi^a™, lorsque m est superieur a l r et comme 
20 etant 1' intervalle A lorsque m vaut 1, 

3) on choisit m valeurs initiales complexes d' un indice 
de refraction complexe n*=n+jk aux m points (Xi, i allant 
de 1 a m, 

4) lorsque m est different de 1, on choisit une loi 
25 d 1 interpolation qui permet de calculer 1 T indice de 

refraction n( a) du materiau sur 1' intervalle B, a 
partir des points (a ir ni) , avec ni=n(ai), i allant de 1 
a m, et lorsque m vaut 1, n(cc) est pris egal au nombre 
ni (ai) sur tout 1' intervalle B, 
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5) on choisit M parametres variables, M etant inferieur 
ou egal a 2m+l, 

6) on choisit une fonction d'erreur Er C¥ , ¥ ) qui 
caracterise l'ecart entre un spectre mesure T et un 

5 spectre experimental Y , 

7) a l f aicle d'une fonction de minimisation de Er(^F r X Y ) 
a M parametres, on effectue la serie d'etapes 
suivante : 

a) a l'aide de la loi d' interpolation des (a^ni) sur 
10 . l'intervalle B, on deduit n(a), a appartenant a B, 

b) & l'aide de n(a) et de l'epaisseur s de la couche, 
et de methodes de calcul de spectres, on calcule un 
spectre theorique In (a) , e) r 

c) on compare Y et Y a l'aide de Er OF, W) et, si 
15 Er C¥, W) est suffisamment petite, c' est-a-dire 

inferieure a une valeur predefinie e, ou est minimale r 
on passe a 1'etape e) , sinon on passe a 1'etape d) , 

d) on fait varier les M parametres variables de fa^on a 
tendre vers le minimum de Er C¥, Y ) et l'on retourne a 

20 1' 6 tape a) , 

e) si Er (Y, W) est inferieure a e, on obtient done un 
ensemble de M parametres variables, pour lequel Er C¥ , 
¥ (n (a,M) ,e) ) est minimum et l'indice de refraction est 
alots pris egal a celui qui a ete obtenu en dernier, et 

25 si Er( x F, W) est superieure ou egale a e on passe a 
1'etape 8) , 

8) on augmente le nombre m de valeurs initiales de la 
fonction a et l'on retourne a 1'etape 2) . 
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2. Procede selon la revendication 1, dans 
lequel on augmente le nombre de valeurs initiales de la 
fonction a en ajoutant une ou plusieurs valeurs aux 
valeurs initiales existantes. 
5 3. Procede selon la revendication 1, dans 

lequel on augmente le nombre de valeurs initiales de la 
fonction a en remplagant les valeurs initiales 
existantes par de nouvelles valeurs initiales dont le 
nombre est superieur au nombre des valeurs initiales 
10 existantes. 

4. Procede selon l'une quelconque des 
revendications 1 a 3 f dans lequel chaque loi 
d r interpolation est choisie parmi les lois 
d f interpolation lineaires, les lois d 1 interpolation 

15 cubiques, les lois d 1 interpolation polynomiales et les 
lois d l interpolation par exemple de type fonction 
spline . 

5. Procede selon l'une quelconque des 
revendications 1 ci 4 , dans lequel les valeurs 

20 initiales de la fonction a sont reparties de maniere 
uniforme dans 1' intervalle A, la repartition des noeuds 
etant ainsi homogene. 

6. Procede selon l'une quelconque des 
revendications 1 a 5, dans lequel a (A) est choisi parmi 

25 X, 1/X et hc/X r ou h est la constante de Planck et c 
est la vitesse de la lumiere dans le vide. 

7. Procede selon l'une quelconque des 
revendications 1 a 6, dans lequel on mesure l'erreur, k 
l'etape 6), sur un intervalle d' interet C qui est 

30 inclus dans 1' intervalle B ou egal a cet intervalle B. 
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8. Procede selon l'une quelconque des 
revendications 1 a 7, dans lequel les M parametres 
variables sont les parties reelles des indices de 
refraction aux points ai, i allant de 1 a m. 
5 9. Procede selon l'une quelconque des 

revendications 1 a 7, dans lequel les M parametres 
variables sont les parties imaginaires des indices de 
refraction aux points a if i allant de 1 am. 

10. Procede selon l'une quelconque des 
10 revendications 1 a 7, dans lequel les M parametres 
variables sont constitues par l'epaisseur du materiau 
dont on cherche l'indice de refraction* 
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